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Greater Impact of Melanocortin-4 Receptor Deficiency 
on Rates of Growth and Risk of Type 2 Diabetes 
During Childhood Compared With Adulthood in 
Pima Indians 

Marie S. Thearle, Yunhua L. Muller, Robert L. Hanson, Meghan Mullins, Maryam AbdusSamad, 
John Tran, William C. Knowler, Clifton Bogardus, Jonathan Krakoff, and Leslie J. Baier 



Features of melanocortin-4 receptor (MC4R) deficiency have 
been observed to be more pronounced in childhood. Longitudinal 
data from a population-based study were used to separate the 
phenotypic effects of MC4R deficiency during childhood and 
adulthood. The MC4R exon was sequenced in 6,760 individuals 
of predominantly Pima Indian heritage, and discovered muta- 
tions were functionally assessed in vitro. Effects on BMI, height, 
and slope of BMI change were assessed during childhood (ages 
5-20 years) and adulthood (ages 20-45 years). Six mutations 
affecting MC4R function, including three that may be private 
to Pima Indians, were found in 159 individuals (2.4%). The slope 
of BMI increase was greater in individuals carrying an MC4R 
mutation compared with noncarriers during childhood but not 
during adulthood. The final adult height obtained was higher in 
individuals with MC4R deficiency. There was an increased risk 
for developing type 2 diabetes in individuals with a defective 
MC4R during childhood and adulthood, but this was only inde- 
pendent of BMI in childhood. The greater rates of body mass 
accumulation and risk of type 2 diabetes before the age of 20 
years in individuals with MC4R deficiency indicate that the 
effects of these mutations are more apparent during the active 
growth of childhood. Diabetes 61:250-257, 2012 




The melanocortin system in humans integrates op- 
posing signals from the hypothalamic pathways 
that regulate food intake and energy homeostasis. 
The major interaction point for the anorexigenic 
signal, melanocyte-stimulating hormone (a-MSH), and the 
orexigenic signal, agouti-related peptide, is the melanocortin-4 
receptor (MC4R). MC4R haploinsufhciency is recognized 
as the most common monogenic form of obesity, with a 
prevalence of 1.7-3.0% in obese individuals (1-6). Children 
carrying an MC4R mutation have been observed to be obese, 
hyperphagic, and hyperinsulinemic, with increased linear 
growth (6-8). Many of these phenotypic features seem to 
become less prominent with progression into adulthood 
(3,7-10). 

Pima Indians living in Arizona have a high prevalence 
of obesity (11) and diabetes (12) attributable to genetic 
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(13,14), environmental, and lifestyle influences (12,15). 
After the first month of life, Pima children have higher 
BMI z scores than the general U.S. population at all ages 
(16). Even during childhood, the diabetes that occurs in 
this population has been shown to consist primarily of 
type 2 diabetes (17). Previous sequencing of MC4R in 
severely obese Pima Indians identified the missense 
change, R165Q, which is known to cause partial loss of 
function (7,18,19), and an A insertion at nucleotide 100 
that causes a frameshift mutation that results in a pre- 
mature STOP at codon 37 (D37stop) (20) prior to the li- 
gand binding and transmembrane domains. To date, this 
latter mutation has not been reported in other populations. 
The prevalence of these two mutations in full-heritage 
Pima Indians is 1.8% (4). Our aim was to extend these 
findings by sequencing the entire coding region of MC4R 
in all participants of our longitudinal survey of health. 
Knowledge of every MC4R coding mutation in the Gila 
River Indian community, without prior selection for obese 
individuals, permitted us to determine the full prevalence 
and to assess the phenotypic impact of MC4R deficiency 
during childhood and adulthood. Based on previous obser- 
vations reported in the literature (6-8), we hypothesized 
that MC4R deficiency may lead to more pronounced effects 
during childhood. 

RESEARCH DESIGN AND METHODS 

Between 1965 and 2007, all members of the Gila River Indian Community aged 
>5 years were invited to participate in the National Institutes of Health lon- 
gitudinal study of health as previously described (11). Individuals participated 
as frequently as every 2 years in examinations that included measurements of 
weight, height, and fasting and 2-h glucose and insulin concentrations during 
a 75-g oral glucose tolerance test (OGTT). The degree of Pima heritage was 
determined by the self-reported number of great grandparents of full Pima 
heritage. The evaluated population was limited to 6,760 people of >50% Pima 
or related Tohono O'odham heritage to reduce confounding from ethnic dif- 
ferences in BMI and risk of type 2 diabetes. Blood pressure measurements 
were available for 5,957 individuals, and self-reported history of menarche was 
available for 3,334 female subjects. All adults provided written informed 
consent; a parent or guardian provided informed consent for children aged 
£18 years, and the child assented. This study was approved by the institu- 
tional review board of the National Institute of Diabetes and Digestive and 
Kidney Diseases. 

Direct sequencing of the MCtyR exon. The entire coding region of MC4R 
was PCR amplified with FideliTaq DNA polymerase (Affymetrix, Santa 
Clara, CA), using primers 5 '-GCAATTTTAGCCTCACAACTT-3 ' (forward) and 
5'AAATCCACAGTGCCTACAACC-3' (reverse). Direct sequencing was con- 
ducted in all subjects using a Big Dye terminator kit (Applied Biosystems, 
Carlsbad, CA) and electrophoresed on an automated DNA capillary se- 
quencer (model 3730; Applied Biosystems). 

Functional characterization of MC4R variants. The open reading frame for 
the common sequence or the variant MC4R DNA was PCR amplified from the 
genomic DNA of previously sequenced subjects and cloned in the vector 
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TABLE 1 



Number of study visits per 5-year age range in control subjects (individuals without an MC4R coding mutation) compared with 
individuals carrying an MC4R loss-of-function mutation 











Age at study visit (years) 








Total 


5-9 


10-14 


15-19 


20-24 


25-29 


30-34 


35-39 


40-45 


Number of study visits 


Control subjects 


3,305 


5,448 


4,390 


3,378 


2,787 


2,400 


1,886 


2,240 


25,834 


MC4R mutation carriers 


80 


110 


101 


72 


64 


42 


30 


52 


551 



pcDNA3.1/V5-His using a TOPO TA expression kit (Invitrogen, Carlsbad, CA). 
Vector with no insert served as a negative control. All constructs were se- 
quenced as described above to confirm the correct variant. One microgram of 
each expression construct was transiently cotransfected with a 250-ng mixture 
of cAMP-responsive element (CRE)-mediated firefly luciferase construct and 
constitutively expressed Renilla luciferase construct (40:1) (SABioscience, 
Frederick, MD) using lipofectamine LTX (Invitrogen) into HEK293 (human 
embryonic kidney) cells maintained in Eagle's minimum essential medium 
supplemented with 10% FBS, 100 units/mL penicillin, and 100 ^g/mL strepto- 
mycin at 37°C with 5% C0 2 . Six hours posttransfection, cells were serum 
starved for 18 h in OPTI-MEM1 plus 0.5% FBS. Graded concentrations of 
a-MSH (Sigma-Aldrich, St. Louis, MO) were then added to the medium for 
a final concentration of 10 -11 to 10~ 6 mol. After 6 h, cells were harvested and 
cAMP-induced luciferase activity, reflecting the level of MC4R activation, was 
assessed using the Dual-Luciferase Reporter Assay System (Promega, Madison, 
WI) on a luminometer (Turner Designs, Sunnyvale, CA). Relative luciferase 
activity was expressed as a ratio of firefly luciferase activity by cAMP to Renilla 
luciferase activity to control for transfection efficiency and normalized to 
a percentage of maximum wild-type activity. Six replicate experiments were 
performed per construct. 

Statistics. The a level for significance was set at 0.05. Comparisons between 
the EC 50 , which is the concentration of a-MSH required to achieve 50% of the 
maximal effect, and maximum activity of the variants were assessed with one- 
way ANOVA followed by Dunnett posttest after nonlinear regression fitting 
of a four-parameter sigmoid dose-response curve using GraphPad Prism 5 
(GraphPad Software, La Jolla, CA). If a variant led to either (i) a decrease in 
the maximum response and/or (Z) an increased EC50, it was considered to be 
a loss-of-function mutation. 

Comparison of rates of growth. For all individuals with two or more 
nondiabetic, nonpregnant examinations during childhood, defined as aged 5-20 
years (n = 3,397; 69 subjects with a loss-of-function MC4R mutation), or 
adulthood (n = 2,645; 52 subjects with a loss-of-function MC4R mutation), the 
slope of BMI change was calculated. The age of adults was restricted to >20 
and :=45 years because the number of exams after the age of 45 years was low. 
To verify that MC4R deficiency was associated with childhood effects in- 
dependent of early-onset excess adiposity, an analysis (n = 276) was con- 
ducted using three control subjects matched by sex, to within 1 year of the 
initial age, and to within 0.1 of the initial BMI z score for each individual with 
MC4R deficiency. Linear regression models were used to compare slopes, with 
sex, degree of Pima heritage, and birth year (to account for secular changes in 
growth) as covariates. 

z Scores and mixed models. Data from individuals without an MC4R mu- 
tation were used to calculate means and SDs for every age by sex for BMI and 
height. Pima population-specific z scores were calculated for each individual 
by subtracting the mean from the individual's value and dividing by the SD. 
For comparison with reference values, BMI 2 scores also were calculated from 
the standards published by the National Center for Health Statistics of the 
Centers for Disease Control and Prevention (21) using the recommended 
equation for normalized z scores. The number of nondiabetic, nonpregnant 
study exams with data for height and weight ranged from 1 to 17 in 5,977 
individuals for a total of 26,385 examinations between 1965 and 2007 (Tables 1 
and 2). Of these, 13,364 examinations occurred in 4,537 children. 

Using SAS software (SAS Institute, Cary, NC), mixed models using a first- 
order autoregressive error structure for repeated measures were used to 
compare individuals with and without a defective MC4R, controlling for degree 
of Pima heritage and birth year as fixed effects and sibling relationship as 
a random effect. Analyses included height and BMI z scores, BMI, final adult 
height, and systolic and diastolic blood pressure. Excluding the z score 
analyses, all models also included age and sex as covariates. Blood pressure 
models were further adjusted for BMI and height. A subset analysis (n - 436) 
for final adult height with three control subjects matched as above by BMI z 
score, sex, and age for each individual with MC4R deficiency was conducted 
to separate the effects of excess adiposity and MC4R deficiency. Analyses 



were adjusted for multiple comparisons using Dunnett posttest whenever 
carriers of MC4R variants were considered as separate groups. Linear regression 
was used to compare the age of menarche between girls with and without MC4R 
deficiency, controlling for birth year and degree of Pima heritage. 
Rates of type 2 diabetes. Time to development of diabetes was calculated 
from the first nondiabetic visit to the onset of diabetes determined from the 
OGTT results (22) or chart review if diabetes was diagnosed between study 
visits. Cox proportional hazards models were used to determine the hazard 
rate ratio (HRR) associated with MC4R deficiency for developing diabetes 
before the age of 20 years with and without BMI as a time-dependent variable 
and including sex, birth year, and degree of Pima heritage as covariates (n = 
4,518; 103 subjects with MC4R deficiency). Similar analyses were conducted to 
assess the HRR for developing type 2 diabetes between the ages of 20 and 45 
years (n = 4,289; 102 subjects with MC4R deficiency). 

RESULTS 

Functional analysis and prevalence of MC4R mutations. 

Sequencing of MC4R in 6,760 (44% male) individuals iden- 
tified one nonsense (D37stop) and nine missense (R165Q, 
R165G, A303P, G231S, I269N, I251L, P275S, V103I, and 
L304F) variants. In vitro analysis revealed that 6 of the 
10 changes caused decreased MC4R function (see Table 
3 and Fig. 1). R165Q, G231S, and I269N were previously 
known to lead to partial loss of function (1,7,23,24), con- 
sistent with our findings. The frameshift mutation, D37stop, 
predicted termination of the protein prior to the ligand 
binding domains yet had some detectable activity with 
exposure to high concentrations of ligand. Two missense 
mutations, R165G and A303P, also had partial activity. In 
contrast, I251L, P275S, V103I, and L304F had in vitro activity 
similar to wild-type MC4R, and individuals with these var- 
iants (n = 39) were excluded from all additional analyses. 

The prevalence of the six loss-of-function MC4R muta- 
tions was 2.4%, with 156 heterozygous individuals and 3 
siblings homozygous for the D37stop mutation. To increase 
the statistical power of additional analyses, all individuals 
with any of the six loss-of-function MC4R mutations were 



TABLE 2 

Number of times control subjects (individuals without an MC4R 
mutation) and MC4R loss-of-function mutation carriers were seen 
for a nondiabetic study visit 



Number of 


Control subjects 


MC4R mutation 


nondiabetic visits 


(n) 


carriers (n) 


1 


973 


27 


2 


918 


21 


3 


817 


20 


4 


723 


17 


5-9 


1,988 


39 


10-14 


384 


12 


>15 


38 


0 


Total 


5,841 





Visits could occur as often as every 2 years. 
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TABLE 3 



Comparison of the in vitro activity (defined by EC50 and maximum response) for each variant identified by sequencing the MC4R exon 





EC 50 (95% CI) 




Maximum response (%) 




Childhood BMI 




Childhood BMI 






Variant 


(nmol/L) 


P 


(95% CI) 


P 


(kg/m 2 ) 


P 


z score 


P 


n 


Control 


o.£ (iL.iL— 4.0) 




1 AA A /'ACC O 1 AO 1\ 

iuu.u (yo.o— lUO. I) 




O/l f -4- A 1 
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<-U.UUl 


Otl O fQI A 1 1 A K\ 

y/.Z (po.4— 114.0J 


A AO 


OQ A -»- A Q 

zy.y _ u.y 


<-U.UUl 


A QA + n 1C 

u.yo _ u.io 


<-U.UUl 


A£ 
DO 


KlODlr 


PCA 7 /1Q ^ nnn CC^ 

oo. I (io.4— Zoy.o) 


^ a An 1 
<U.UU1 


QQ 1 /'AT A 1 AO 

oo. 1 (0/.4— lUz.o) 


A in 
U. Id 


OA "7^11 

zy. ( ± i.i 


^--A AA1 
<U.UU1 


A Q7 -t- A 1 O 


-~-A AA1 

<U.UU1 


do 


A303P 


10.8 (3.9-33.5) 


0.05 


51.9 (45.7-62.1) 


<0.001 


32.5 ± 2.0 


<0.001 


1.16 ± 0.33 


<0.001 


10 


D37stop 


104.2 (44.1-246.6) 


<0.001 


52.4 (42.5-66.4) 


<0.001 


31.2 ± 1.2 


<0.001 


1.22 ± 0.21 


<0.001 


30 


G231S 


9.3 (5.5-15.6) 


0.006 


110.3 (100.8-123.8) 


0.09 


25.7 ± 2.0 


0.6 


0.28 ± 0.33 


0.4 


12 


I269N 


9.8 (5.7-16.8) 


0.01 


103.8 (93.0-118.6) 


0.35 


21.9 ± 3.2 


0.4 


-0.44 ± 0.54 


0.4 


7 


125 1L 


3.6 (2.2-6.0) 


0.71 


100.3 (94.7-109.9) 


0.53 


20.0 ± 1.7 


0.007 


-0.74 ± 0.29 


0.01 


14 


P275S 


4.7 (2.7-8.3) 


0.56 


99.4 (92.6-110) 


0.74 


22.2 ± 5.9 


0.7 


-0.29 ± 0.97 


0.8 


2 


V103I 


5.2 (3.1-8.9) 


0.31 


95.2 (90.1-104.3) 


0.44 


24.7 ± 1.5 


0.9 


-0.01 ± 0.24 


0.9 


18 


L304F 


6.6 (3.5-12.5) 


0.12 


98.3 (91.5-109.2) 


0.94 


23.5 ± 3.1 


0.9 


-0.12 ± 0.64 


0.9 


5 



For each variant, the mean childhood BMI and BMI z score is given for all individuals from the study population found to carry the respective 
MC4R variant. The variants are listed in the same order as portrayed in Fig. 1. 



combined into one group (individual results for the six 
single nucleotide polymorphisms [SNPs] are presented in 
the Supplementary Table). By the age of 20 years, 81% of 
children with an MC4R loss-of-function mutation were 
considered to be at risk for obesity, denned by a BMI 
greater than the U.S. population's 95th percentile, com- 
pared with 50% of control subjects (x 2 = 41; P < 0.0001). Of 
children at risk for obesity, 3.6% had MC4R deficiency. By 
the age of 45 years, 90% of individuals with MC4R de- 
ficiency were obese, denned by a BMI >30 kg/m 2 , versus 
76% of control subjects (x 2 = 9; P < 0.002). Of obese 
adults, 2.8% carried an MC4R loss-of-function mutation. 
Differences in rates of BMI gain. Among the 4,537 indi- 
viduals with at least one measurement for BMI in childhood, 
106 had a loss-of-function MC4R mutation (40 with R165Q, 
26 with R165G, 21 with D37stop, 8 with A303P, 8 with 
G231S, and 3 with I269N). The average age of initial child- 
hood evaluation was similar for individuals with and with- 
out an MC4R mutation (aged 10.9 ± 3.7 vs. 11.1 ± 4.2 years; 
P = 0.6). The ratio of the sexes and the degree of Pima 
heritage was similar between the two groups. In children 
with more than one evaluation (n = 3,397), the rate of 



growth, as measured by the slope of BMI increase, was 
1.50 ± 0.08 kg/m 2 /year in children heterozygous for a loss- 
of-function mutation compared with 1.18 ± 0.01 kg/m 2 / 
year in control subjects (P < 0.001) (Fig. 2). Neither the 
average number of observations (3.7 ± 1.4 vs. 3.8 ± 1.4; 
P = 0.6) nor the age range between the first and last ob- 
servation (aged 8.0 ± 3.4 vs. 8.1 ± 3.6 years; P = 0.9) dif- 
fered between the two groups. In subset analysis with 
control subjects matched by BMI z score, age of initial 
visit, and sex, the degree of difference between the groups 
was attenuated but still significant (1.55 ± 0.08 vs. 1.35 ± 
0.01 kg/m 2 /year; 95% CI for the difference 0.03-0.37; P = 
0.02). Children with a loss-of-function MC4R mutation had 
a higher adjusted BMI z score over all childhood exami- 
nations (0.89 ± 0.10 vs. 0.00 ± 0.02; 0.70-1.09; P < 0.001). 
Corresponding BMIs were 29.7 ± 0.6 kg/m 2 and 24.9 ±0.1 
kg/m 2 for those with and without MC4R deficiency, re- 
spectively (3.5-5.9; P < 0.001) (Fig. 3). For comparison, 
U.S. population BMI z scores would be 2.01 ± 0.78 for 
children with MC4R deficiency and 1.32 ± 0.99 for the rest 
of the population. The three siblings homozygous for 
D37stop had an average BMI of 47.2 ± 3.4 kg/m 2 (40.5-53.8; 



A B C 




Log a-MSH (M) Log a-MSH (M) Log a-MSH (M) 

FIG. 1. Functional characterization of the identified MC4R missense/nonsense variants. MC4R activity (means ± SE) in response to a logarithmic 
increase of a-MSH is expressed as a percentage of the maximum wild-type MC4R activity. The best-fit estimate of the EC 50 , the concentration of 
a-MSH required to achieve 50% of the maximal effect starting from basal activity, and the 95% CIs were obtained from nonlinear regression fitting 
of the sigmoid dose-response curve using GraphPad Prism 5. A functional MC4R defect was defined as either a statistically significant decrease in the 
maximum response or a significant increase in the EC 50 . Those variants showing the greatest in vitro effects are shown in A, variants showing 
a milder effect are shown in B, and those variants with no detectable difference from wild-type MC4R are shown in C. Vector pcDNA3.1 served as 
a negative control for background cAMP accumulation. 



252 DIABETES, VOL. 61, JANUARY 2012 



diabetes.diabetesjournals.org 



M.S. THEARLE AND ASSOCIATES 



p<0.001 




Childhood 


Adulthood 


p<0.001 








Childhood 



Adulthood 



FIG. 2. A: Rates of BMI change in kilograms per meters squared per 
year (kg/m 2 /year) in individuals without diabetes with a loss-of-function 
MC4R mutation (■) and with normal function (□) during childhood 
(ages 5-20 years; n = 3,397, 69 subjects with MC4R deficiency) and 
adulthood (ages 20-45 years; n = 2,695, 52 subjects with MC4R de- 
ficiency), adjusted for birth year, degree of Pima heritage, and sex. B: 
Rates of BMI increase in kilograms per meters squared per year (kg/m 2 / 
year) in childhood and adulthood in 1,501 individuals who were present 
in both the adult and childhood datasets adjusted for birth year, degree 
of Pima heritage, and sex. Twenty of these individuals had MC4R de- 
ficiency (H). 



P < 0.001) and a BMI z score of 4.5 ± 0.56. As shown in 
Fig. 3, the average BMI at almost every age for girls and 
boys with MC4R deficiency was greater than the 50th 
percentile on the Pima population-specific growth chart 
and greater than the 95th percentile on the U.S. population 
growth chart. 

The degree of difference in the MC4R loss-of-function 
group stabilized after the age of 20 years, such that the 
adjusted z score for adults was 0.59 ± 0.11 in people with 
MC4R deficiency (95% CI for the difference 0.31-0.74; P < 
0.001), corresponding to an average BMI of 38.9 ± 0.9 kg/m 2 
compared with 34.7 ±0.1 kg/m 2 in control subjects (2.5-5.9; 
P < 0.001). The slope of the BMI increase was no longer 
greater in adults with a defective MC4R (0.47 ± 0.08 vs. 
0.43 ± 0.01 kg/m 2 /year; P = 0.7) (Fig. 2). We performed 
subset analyses to verify the differences observed between 



adults and children. These analyses, which yielded simi- 
lar results, were (1) to restrict the population to full- 
heritage Pima Indians (n = 3,188), (2) to exclude children 
who developed type 2 diabetes (n - 5,838), (3) to include 
only individuals with multiple visits during both child- 
hood and adulthood (n = 1,501; Fig. 2E), (4) to limit the 
analysis to individuals with three or more visits (n = 2,650), 
and (5) to restrict childhood age to between 5 and 16 years 
(n = 2,877). 

Risk of type 2 diabetes. Children carrying an MC4R mu- 
tation had higher fasting glucose (5.09 ± 0.04 vs. 4.95 ± 0.01 
mmol/L; P < 0.001), 2-h glucose (5.81 ± 0.09 vs. 5.56 ± 0.02 
mmol/L; P = 0.009), fasting insulin (279.6 ± 15.6 vs. 
231.6 ± 2.4 pmol/L; P < 0.001), and 2-h insulin concen- 
trations (1,129.2 ± 64.8 vs. 962.4 ± 9.6 pmol/L; P < 0.001); 
however, after accounting for BMI, none of these differences 
remained statistically significant. Likewise, differences be- 
tween nondiabetic adults with or without a defective MC4R 
in fasting or 2-h glucose or insulin concentrations were not 
significant after adjustment for BMI. Fifteen individuals 
(10.1%) with MC4R deficiency developed type 2 diabetes 
before the age of 20 years compared with 154 (2.6%) 
control subjects. The HRR, adjusted for sex, birth year, 
and the percentage of Pima heritage, for developing di- 
abetes before the age of 20 years in individuals carrying 
a loss-of-function MC4R mutation was 4.3 (95% CI 2.3-8.3; 
P < 0.001) (Fig. 4). The HRR was attenuated but still 
significant after further adjustment for BMI (3.3 [1.2-9.2]; 
P = 0.03). MC4R deficiency also increased the risk of de- 
veloping diabetes between the ages of 20 and 45 years, with 
an HRR of 2.0 (1.4-2.7; P < 0.001) (Fig. 4). However, after 
adjustment for BMI, this risk was fully attenuated (1.4 [0.9- 
2.2]; P = 0.12). 

Height differences. Final adult height was higher in the 
group with MC4R deficiency (168.8 ± 0.5 vs. 167.2 ± 0.1 
cm [95% CI for the difference 0.6-2.6]; P = 0.001). In ad- 
dition, the average adjusted height z score for those with 
MC4R deficiency was significantly higher during childhood 
at 0.25 ± 0.09 (P = 0.01) (Fig. 5). Differences between 
groups were similar in the subset analysis, with control 
subjects matched by BMI z score, age, and sex. The three 
siblings homozygous for the D37stop mutation had a higher 
average adult height of 175.0 ± 2.9 cm. 
Blood pressure differences. Before the age of 20 years, 
MC4R deficiency was associated with decreased systolic 
blood pressure in the loss-of-function group even after 
adjusting for BMI, height, sex, age, and degree of Pima 
heritage (106.3 ± 0.8 vs. 108.8 ± 0.1 mmHg [95% CI for 
the difference 0.8-4.1]; P = 0.003) but not with a differ- 
ence in diastolic blood pressure (P = 0.7). No differences 
were observed between groups in either systolic (P = 0.4) 
or diastolic (P = 0.2) blood pressure during adulthood 
after adjustment for the above confounders. 
Differences in menarche. In the 3,334 women (79 with 
MC4R deficiency) with data recorded, the presence of an 
MC4R mutation was associated with a later onset of 
menses by ~6 months (12.5 ± 0.2 vs. 12.0 ± 0.03 years [95% 
CI for the difference 0.12-0.83]; P = 0.008), even after 
adjusting for birth year and degree of Pima heritage. In- 
clusion of BMI in the model did not alter the results. 



DISCUSSION 

This study describes a detailed investigation of the prev- 
alence of MC4R coding mutations, assessment of their 
functionality, and analysis of fife course growth patterns in 
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FIG. 3. A: Pima population-specific BMI z scores in children (age 5-20 years) with (■) and without (□) partial loss-of-function MC4R mutations 
after adjustment for birth year, sibling relationships, and degree of Pima heritage. B: Mean BMIs of children with (H) and without (□) loss-of- 
function MC4R mutations. BMI analysis included age, sex, sibling relationships, degree of Pima heritage, and birth year as covariates. C: The BMI 
growth chart for girls shows that average BMI for individuals with MC4R deficiency (#) was consistently greater than the 50th percentile for the 
Pima population (gray diamonds, dashed lines) and greater than the 95th percentile for the U.S. population (white boxes, dotted lines) at every age. 
D: Corresponding boys' growth chart for BMI. Gray diamonds, dashed lines = 95th, 50th, and 5th percentile for the Pima population; white boxes, dotted 
lines = 95th and 50th percentile for the U.S. population. (A high-quality color representation of this figure is available in the online issue.) 



a longitudinal study of health with subjects of predomi- 
nantly Pima Indian heritage. In this founder population, 
10 missense/nonsense mutations were found, of which 6 
cause loss of function in vitro. Three of these mutations, 
D37stop, R165G, and A303P, have not been described in 
others and may be private to this population. Although 
each mutation individually was rare, 2.4% of the population 
was heterozygous for MC4R deficiency, providing support 
for the hypothesis that accumulation of multiple rare 
mutations in a population may have a role in susceptibility 
for a common disease. It is well established that mutations 
in MC4R that result in decreased cAMP production lead to 
an increased risk for obesity in both rodents and humans 
(2,7-10,19), which seems to be attributed to hyperphagia 
as well as to decreased energy expenditure (4), likely as 
a result of a reduction in sympathetic nervous system 
activity (25). Our results show that the phenotypic effects 
of MC4R deficiency are more pronounced during child- 
hood. 

Previous investigators have anecdotally noted that the 
phenotype of people with MC4R deficiency may differ be- 
tween children and adults (5,8). For example, the hyper- 
phagia seen in childhood has been reported to be less 
pronounced in adults, abating sometime in the late teenage 
years (5,7,8). Our work adds to this previous literature by 



using data from a large, longitudinal population study to 
confirm in a more statistically robust way the differential 
effects of a defective MC4R. Pima children carrying a loss- 
of-function MC4R mutation had greater rates of weight 
gain, leading to higher BMIs. After completion of child- 
hood growth, rates of BMI gain stabilized and were no 
longer different from the rest of the population. Thus, the 
BMI difference between those with and without MC4R 
deficiency did not increase further during adulthood but 
stayed stable at ~5 kg/m 2 . As this was a population-based 
study, there was no selection bias for obese individuals. 
There is, however, a high prevalence of obesity in the Pima 
population during both childhood and adulthood 
(11,15,26). Previous reports conflict on the effect of car- 
rying a defective MC4R on the rate of BMI gain in adults, 
with one study (10) reporting no increased BMI gain per 
year and another (6) demonstrating an increased pene- 
trance of obesity at age 40 years compared with self- 
reported BMI at age 20 years. However, no other study 
has reported prospective standardized weight and height 
measurements to investigate separate childhood versus 
adult effects of MC4R deficiency. 

The reasons for accelerated weight gain only during 
childhood, and not adulthood, are not clear. Hyper- 
responsiveness of the pathways that frame the appetitive 
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FIG. 4. Survival curve for time to type 2 diabetes between the ages of 
5 and 45 years. By age 20 years, 2.6% of the population without an 
MC4R mutation (solid line) developed type 2 diabetes compared with 
10.1% of those with MC4R deficiency (dashed line). The increased risk 
for type 2 diabetes in those with MC4R deficiency continued into 
adulthood, although it was less pronounced. The Cox proportional 
hazards analysis included the following covariates: degree of Pima 
heritage, birth year, and sex (P < 0.001). 

response to low peripheral energy stores during childhood 
may have had greater survival benefit in times of food 
deprivation. There also may be pathways that interact with 
the melanocortin system that are less active after the 
completion of linear growth. In support of the first pos- 
sibility, a recent study showed that carrying a greater 
number of SNPs associated with adult BMI decreased the 
risk of failure to thrive during early infancy (27). In 
support of the latter possibility, MC4R has been shown to 
be important in the regulation of the growth hormone 
axis (28), which is known to decrease in physiologic 
importance with progression to adulthood. A common SNP 
near MC4R also has been shown to have a stronger asso- 
ciation with BMI before the age of 20 years, which weak- 
ened during adulthood (29). 

MC4R deficiency has been reported to lead to an early 
increase in stature (7,8,28), although not every study has 
confirmed this finding (2,30). Our study found an increase 
in childhood height z score in individuals with MC4R de- 
ficiency as well as a taller final height by -1.5 cm. The 
effect on linear growth may be related to the severity of 
MC4R deficiency, which would explain why our study, with 
partial loss-of-function mutations, found a relatively small 
effect and another recent study demonstrated a larger effect 
on final adult height of 5-7 cm in individuals heterozygous 
for complete loss-of-function mutations (28). Pima Indian 
children have a different pattern of linear growth than the 
general U.S. population and exhibit increased early growth 
such that, on average, they are taller than the overall U.S. 
population between the ages of 4 and 14 years. However, 
final adult height in the Pima population is similar to the 
general population (16). 

Individuals with the I251L genotype had a significantly 
lower average BMI than the group with a normal MC4R 
genotype. We did not find any difference in in vitro signaling 
between this variant and the wild-type MC4R. Previous 
studies have found that the I251L polymorphism may be 
protective against obesity (31). Our association results seem 
to indicate a protective effect as well; however, it should be 



noted that the group of individuals with this variant included 
no full-heritage Pima Indians and had an overall lower per- 
centage of Pima heritage than the general population, which 
also may account, in part, for their lower mean BMI. 

The lower blood pressure reported by others in adults 
(25) was only observed during childhood in our study. Our 
exclusion of visits after the development of type 2 diabetes 
or the inclusion of all subjects regardless of obesity status 
may account for the difference. Our childhood results, as 
well as the results of others (25,32), indicate that MC4R 
deficiency dissociates the link connecting hypertension 
with obesity and raises the possibility that the increased 
sympathetic nervous system activity seen with excess ad- 
iposity may be mediated through MC4R. The delayed men- 
arche observed in this study may be related to a similar 
phenomenon, as the sympathetic nervous system is impli- 
cated in increasing gonadotropin-releasing hormone secre- 
tion prior to pubertal onset, and women with familial 
dysautonomia also have been reported to have delayed 
menarche (33). The delayed menarche was an unexpected 
finding because it is contrary to what is usually observed 
with increased adiposity; however, it is interesting that, 
whereas other common obesity variants associated with 
early menarche in a genome-wide association study meta- 
analysis, a SNP near MC4R did not (34). 

Consistent with higher BMIs, children with MC4R de- 
ficiency had higher glucose and insulin concentrations 
during an OGTT. However, these differences did not per- 
sist after adjustment for BMI. This is in contrast to the 
classic description of MC4R deficiency, which describes 
hyperinsulinemia out of proportion to BMI that improves 
with age (7,8), although this has not always been con- 
firmed by others (9). We found an increased risk for type 2 
diabetes with MC4R deficiency. Before the age of 20 years, 
this risk was only partially attributable to the increased 
BMI at a young age. Others have not found an increased 
prevalence of type 2 diabetes among adults with MC4R 
deficiency (5,9,25), but we did find a modest risk of de- 
veloping type 2 diabetes during adulthood that was fully 
attributable to the increased BMI. The increased risk for 
type 2 diabetes independent of BMI during childhood was 
possibly attributable, in part, to the rapidity of the weight 
gain. Potentially, the elevated risk is only evident in pop- 
ulations with high susceptibility for developing type 2 di- 
abetes. It is unclear why MC4R deficiency would impart 
additional risk for diabetes beyond that conferred by an 
increased BMI, although recent studies have implicated 
MC4R in potentiating downstream insulin signaling (35,36). 
It also is possible that BMI did not fully represent the effect 
of excess adiposity in developing insulin resistance and type 
2 diabetes in the model. 

In conclusion, in this population-based study, the impact 
of MC4R deficiency on the rate of body mass accumulation 
was greatest during childhood and became similar to the 
rest of the population after completion of childhood growth. 
We also found an increased risk for type 2 diabetes in 
individuals with an MC4R defect, which was more pro- 
nounced during childhood and independent from BMI 
during childhood. In addition, carrying a partial loss-of- 
function MC4R mutation was associated with a taller 
adult height, lower childhood systolic blood pressure, 
and later menarche. Understanding how monogenic changes 
cause differing phenotypic expression during childhood and 
adulthood may help to further insights into common obesity 
and emphasize the need for prevention efforts during 
childhood. 
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GROWTH IN PIMA INDIANS WITH MC4R DEFICIENCY 
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FIG. 5. 4: Final adult height In individuals with a loss-of-function MC4R mutation (» compared with those without an MC4R mutation (□), 
adjusted for sex, birth year, degree of Pima heritage, and sibling relationships. B: Average height z scores during childhood (ages 5-20 years) for 
individuals with a loss-of-function MC4R mutation (■) and for those without an MC4R mutation (□), adjusted for birth year, degree of Pima 
heritage, and sibling relationships. C: The height growth chart for girls shows the average height for individuals with MC4R deficiency at every age 
(#) in comparison with the remainder of the Pima population (gray diamonds, dashed lines). This growth chart also shows the early increased 
linear growth seen in the Pima population compared with the U.S. population (white boxes, dotted lines). D: Corresponding boys' growth chart for 
height. Gray diamonds, dashed lines = 95th, 50th, and 5th percentile for the Pima population; white boxes, dotted linos = 95th, 50th, and 5th 
percentile for the U.S. population. (A high-quality color representation of this figure is available in the online issue.) 
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